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Abstract 

The top quark decay process t — > c T is studied in the Standard Model and the 
Minimal Supersymmetric Standard Model with conserved and violated i?— parity. 
In dealing with the heavy quarkonium, T, production, we take both color-singlet 
and -octet mechanisms into consideration. With numerical analysis it is found 
that within the Standard Model the decay rate of the concerned process is beyond 
the border of reachable level of the next round of experiment, though it is much 
larger than that of the rare top decays, while extra new contributions from beyond 
the Standard Model could enhance the decay rate by several orders. It is also 
noticed that a precise measurement of the charged Higgs mass without tan (3 
dependence is possible via this process. 
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1 Introduction 

The discovery of the top quark in 1995 fl], H §| marked the triumph of the Standard 
Model(SM) of particle physics. Since the discovery, more interest on researches in top 
quark physics is stimulated in both experimental and theoretical respects. Because of 
its large mass which is close to the electroweak scale, the top quark physics research 
may play an important role in the study of the electroweak symmetry breaking and 
therefore of the origin of the fermion masses. Through the study on it, especially its 
rare decays, people also hope to find clues of possible "new physics". 

Within the framework of SM the dominant top quark decay mode is t —>■ bW, which 
has a time scale tw much shorter than that of the non-perturbative QCD effect. This 
character of top quark makes it behaves almost like a free particle, which is helpful to 
the precise study of it. With the operation of CERN Large Hadron Collider (LHC) 
in the future there will be enormous top quark events obtained every year and that 
provides it practically possible to make a more detailed study on many properties of 
top quark. 

Up to now, the rare top decays have been investigated by several groups within the 
SM@, I and beyond, e.g. in 2-Higgs Doublet Model(2HDM) [g | and Supersymmet- 
ric(SUSY) models In SM, the branching ratios are heavily suppressed because of 
the Glashow-Iliopoulos-Maiani(GIM), typically B{t ^ cH) < 10'^^ §. Though they 
may be greatly enhanced by new contributions from beyond the SM, the rare top decays 
are still quite small with only few cases up to the measurable levels of the near future 
top producing machines, like the upgraded Tevatron and LHC The decay channel 
i — »• cT was once investigated for the aim of studying CP- violation and it was also 
claimed having an accessible rate for the future experiments based on a quite rough 
assumption on the non-perturbative sector of the heavy quarkonium 0, which leads to 
about two orders of magnitude overestimate of the rate, to be seen in the following. 

In this paper, we will investigate the top decay t cT process with taking into ac- 
count both the color-singlet and -octet mechanisms in the T production calculation. We 
find the process is worth being carefully investigated in testing the SM and especially in 
searching new physics signatures for the following arguments. First, comparing with the 
meson decays, the top quark decays possess less theoretical uncertainties, e.g. without 
non-perturbative long distance influence because of its large mass scale. Second, the 
concerned decay mode occurs at tree-level, which may cause the decay rate larger than 
that of rare top decays. Third, theoretically the heavy quarkonium, T, can be treated 
in a relatively high precision by virtue of the Non- relativist ic QCD(NRQCD) factoriza- 
tion scheme |]TU[. Besides, from the experimental point of view, heavy quarkonia have 



clean signatures in their leptonic decays. 

The paper is organized as follows. In section 2 we present the formalism of the 
concerned process t — cT in a general framework. In section 3 the obtained formal- 
ism is applied to some specific models, the SM and Minimal SUSY Standard Model 
(MSSM) with and without i?— parity violation (MSSM=pi?p). The numerical evaluation 
is proceeded in section 4, and in the last section some conclusions are made. 



2 Formalism 



To lowest order, it is quite clear that the interested decay channel can be induced by 
intermediate bosons, for example the W in the SM. In order to be more general and 
convenient to incorporate contributions from different models, in the following we give 
out the expression of the decay amplitude in a universal form. 

The tree-level amplitude of t — > c T could be expressed as 



M = 2V2Gp ml, iV:,Vtb) Y^lr + R 



X 



i=± 



X 



where X stands for the intermediate particles (W^— boson, charged Higgs, etc.), which 
are normalized by the standard VT-boson coupling for the sake of convenience; R repre- 
sents the remnant in the final states in respect to the T production. For abbreviation, 
the caret means the corresponding quantity is normalized by the top quark mass m^, 
i.e. miy = my/ jvnt. The four-quark non-local operators appearing in Eq. (|l]) are of the 
types 



O 



X 



b (Em 



(2) 



(3) 



where g± and h± represent the couplings of the fermion interactions, the Gi and Hi are 
the possible tensor structures of the currents. It is noticed that although in principle 
the quarkonium can be generated through the neutral current of form (|^), however, it 
is obvious that this kind of possible operators, which produce T from a neutral heavy 
vector boson, are of kinematically disfavored ones in leading order since the possible 
large virtuality of intermediate particles. Therefore, it is proper we restrict ourselves 
to the dominant contribution, i.e. form (|^). As well, in our discussion the penguin 
process of the quarkonium production via gluon fragmentation is also dropped because 
that is suppressed by both the GIM mechanism and the smallness of the higher order 



quarkonium Fock states [11 



In the top quark rest frame, the Fix, which denotes the propagator of X, can be 
expressed as Fix = Fi{l,mc,mr,mx,Tx) with 



F+{l,x,y,z,w) = 2[\^x^ -2z^ ^2izw 

F_{l,x,y, z,w) = [y'^ — z^ + izw 



(4) 
(5) 



This is obtained by performing the integrations over the final state phase space at the 
level of decay width, while adding the unstable particle decay width to the propagator, 
and, here and in the after, the lowest order approximation relations in the quarkonium 
system, E-^ = 2Eb and mr = 2mb are taken, or in other words, in the effective La- 
grangian method the relative movement effects of heavy quarks within the quarkonium 



appears only in the nonperturbative sector. Since we only consider the in further 
discussion, we can simply drop the subscript i in Fix, and then Fx will represent F+x 
for short. 

With performing the Fierz transformation on (Q) and considering of only the suitable 
operators in producing T, the operator O'^ could be read as 



where iV„ is the number of colors and 



O 



X 



Y^gfgf [ci,M h^P^h 



i=± 



O 



X 



Y^gtg'r [cl,P^\''t] [h^Piy^^h 
i=± 



(6) 

(7) 
(8) 



with A*^ representing the S'[/(3) Gell-Mann matrices. Here P± = (1 ± 75)/2 are chirality 
projection operators. At the first sight Eq.(^ appears to be the same structure with (|), 
actually they possess different natures in describing quarkonium production. Therefore, 
the operators in Eq.([l|) may be rearranged as 



i=± 



X 



(9) 



where the first term is called the color suppressed one. Next, by taking only the leading- 
order(LO) electroweak effects into account one could replace C± by 



1 , C_ = . 



(10) 



The higher order corrections to these values, either from the QCD or higher order 
electroweak corrections, should be small because of the large top mass scale. Whereas, 
together with the arguments in before, the C_ term can definitely be discarded without 
losing accuracy. 

In the NRQCD framework of factorization, the heavy quarkonium production can be 
properly separated into two sectors, the perturbative QCD calculable part and the non- 
perturbative universal matrix elements. The latter is ordered in the relative velocity, 
f , of the quarks within the quarkonium system. According to the NRQCD, a generic 
5— wave T state can be schematically described by the Fock state decomposition fir 



|T) = 0{l) hh[^S'^']) + 0{v) hh[^Pf^]g) + 0{v 



ric<(8) 



+0{v') gg 



0{v') bbfD^}^'^]gg) + - 



'111 



where the usual spectroscopic notation '^^^^L j is used, and the color states are labeled 
by (1) and (8) superscripts for color-singlet and -octet, respectively. It should be noted 
that in order to express the production rate in terms of NRQCD matrix elements rather 
than the quarkonioum wavefunction, the h and h must be created with a separation small 



compared to the size of the wavefunction, which is of order l/rribV. This condition can 
be expressed as 

as is shown in Fig. 0. Taking the typical value of v"^ ~ 0.08, the value of the RHS of 
eq. ([l2|) is about 0.5 GeV"^ and it is obvious, from the figure, that the entire region of 
mx is allowed for the presumed magnitude of A"s decay width, i.e. ~ 0(l)GeV. 

Based upon preceding preparations the t cT amplitude can now be expressed as 

X |i \c^,p,\-t\ ^ ]^\ci,PA (13) 

with Tl and 7^" stand for the color-singlet and -octet non-perturbative matrix elements. 



q-fl 



T + X 



T + X 



(14) 
(15) 



Here, X represents the possible soft gluon radiations in the process of forming the 
Onium. 

According to expansion (|ll]), the contribution of higher order Fock states evolving 
into the quarkonium by soft gluon(s) emission or absorption are superficially suppressed 
by some orders of the small quantity of v"^ relative to the leading singlet one, however, 
they may play an important role in special cases [Q, where the higher Fock state 
involved processes are compensated by the enhancements originated from other kine- 
matical variables or coupling constant. Notice that the quarkonium production process 
interested in here may have a factorized form provided by the NRQCD, in the following 
for the sake of later convenience the matrix elements in eqs. (pjif) and ([T5|) are squared 
with the soft radiations being summed over and then parameterized as 



X 



(16) 



^(Tr)*T3'^^ = 5'^n(/f^")V^P^ + m|(/f^^)2e;^^e^ + (/f ^i)^ e^-^ 6^1 . (17) 



X 



l;8,2S+iL^ 



Here, the constants ' ' , the "decay constants", stand for the non-perturbative 
parts in the quarkonium producton with a mass dimension, which can be related to the 
vacuum-to- vacuum matrix elements of the NRQCD 4-fermion operators as shown later; 
the pr and are the momentum and polarization vector of T. 

As compared with the color-octet matrix elements, which may now only be deter- 
mined through experiment, the color-singlet part can be determined by several means 



with more accuracy, e.g. from the quarkonium leptonic decay T 
the next to leading order calculation, the decay width is |1T3| 



r(T 



6 nif 




', etc. Up to 



lQas{mr) \ 



where A(l, x) = A(l, x, x) is the triangle function and A(l, x, y) = l+x'^+y'^ — 2{x'^+y'^ + 
x^y^). With experimental values ||I4[ as input, using the running couplings a{mr) = 
1/133, as{iTLr) = 0.178 and mr = 9460.37 ± 0.21 (MeV), one immediately obtains the 
average value of 



/t 



GeV^ 



(19) 



Here, the uncertainty mainly comes from Fcxp. 



The color-octet decay constants defined in Eq. { \[7\j are correlated with the NRQCD 
matrix elements, which are phenomenologically determined from experiments, whereas 
some quantitative diversities of different fit still exist among them. Taking values fitted 
in Ref. |[15||, we have 



if/ 







CxpSf^] 0)/(12mx) = (5.2 ±1.4) x 10-^GeV' 



;o 

2/0 



{8)1 



Oy4mx = (5.2 ± 1.4) x 10"^ GeV^ , 
0)/3m^ = (3.2 ± 2.8) x 10"^ GeV^ . 



(20) 

(21) 
(22) 



Here, (0 



(9x[-^*-^^] 0) are the non-perturbative color-octet matrix elements defined in 



101 . In evaluating the /-J s the approximation of heavy quark spin symmetry has been 



used, i.e. 



(8)1 
J I 



3(0 



10(8)1 



(2J+1)(0 
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(23) 
(24) 



It should be noted that the above calculations on the "decay constants" is just an 
order estimation, especially the color-octet part where large uncertainties still survive. 
Therefore, the results derived in this paper are correspondingly accurate up to about 
an order. 



3 Applying to Models 



With the formalism obtained in preceding section, it is straightforward to make an 
estimation of the branch ratios of the process t — > cT in some specific models. For 
top quark decays, because the overwhelming mode is T(t bW), at least in standard 
model, the branching ratio B{t ^ cT) can be expressed as 



B{t^cT) 



r(t^cT) 



T{t bW) ' 

with the assumption that B{t ^ bW) ^ 1 in all cases in the discussion. 
To next-to- leading order in strong coupling constant, the decay width of t 



16 



T{t bW) 



8^271 



w 



1 - 



2as{mt) /27r^ 
37r V 3 
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(25) 



61^ is 



(26) 



where T'w = J-''{l,mf„rhw) and 

r (1, X, y) = {l- x^f + (1 + x^)y' - 2y^ 



(27) 



Based upon Eq. (|T3|) and considering of the possible color-octet contributions, the 
decay width of t — c T has the form 



57r 



/t 



Ai + 8^(/«)2A, 



X8 



Here, A,- contain contributions from various A"— mediated diagrams, and 



with 



A" 



AT 



(2J 



(29) 

(30) 
(31) 



(ai,a8p5i],a8[i5o]5«8[3Pi]) = (-12,-12,4,-12). 
Here, the auxiliary functions T' j are 

T'\ = J^'spsi] = -^'sppi] = T'{l,mc,mr) , (32) 
which in principle could be complex numbers, and 

rs[^So] = a-^lf-a+^'c)^r- (33) 
In the following, the explicit form of the function A, in concrete models will be given. 



3.1 SM 

In the SM, to leading order the concerned decay process is generated by the FT— boson 
mediated diagram via interaction 



9 



Cw± = ^ V,, [u, 7^P- dj] W+ + h.c. 



(34) 



where Vij denotes the CKM matrix element. The function Aj defined in the previous 
section thereof is 



, SM 



with 



Gj- — , G_ — Fvu . 



(35) 



(36) 



Here, the definitions of and J^'j coincide with that in above discussions. 



3.2 MSSM±i?p 



Now, we try to apply the obtained formalism to models beyond the SM, the MSSM±i?p, 
to see whether there would be some obviously different behaviors. For convenience to 
see the new contributions, i.e. the deviations from aJ^^'^ in the SM, we define 



A. 



Af (1 + 5j) 



where. 



(37) 



2 Re 
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2 ^ ^ 2 

anrricm^ ^ 



1 + 



\ 5G. 



G 



qWs\ 



(38) 



Here in the expression, 5G± contain contributions from A"— mediated diagrams except 
for X = WsM- 



In the MSSM+_Rp[|T^, new contributions at tree- level are induced only by charged 
Higgs bosons through the effective Lagrangian 



C 



9 



Ui 



cot [3 niui P- + t^'^ P+ ) dj + h.c. 



(39) 



where tan /5 = t>2/fi is the ratio of the vacuum expectation values of two Higgs doublets. 
Since the Higgs bosons couple to quarks via the scalar interaction, after the Fierz 
transformation, although in principle there are many kinds of operators, only the mixing 
terms, P±Pzp, would exist for the Upsion production, which consequently leads to a 
cancellation of the tan (3 in couplings, like 



6G^ 



1 rhcThf 



m 



SG_ = 1 ^F^. . 



w 



m 



(40) 



w 



Although the mixing terms are suppressed relative to the squared terms, P^, by a factor 
of m-^/m^r, the nature of independence of the tan/? is of an advantage of it, which may 
lead to a precise measurement of the charged Higgs mass in the future experiments. 
Explicitly, the 5j reads as 



Re 



1 

+4 



1 mx Fh± 
4 FwsM 
1 mx Fh± 



1 + 



4 rhlr Fi 



which behaves like an expansion of 

1 mx Fh± 



4mL F, 



' 9 ao-m^mx 



(41) 



7%x 



(42) 



which tends to be smaller with the increase of the charged Higgs mass in the region 
outside the threshold as shown in Fig. Furthermore, the coefficients contain small 
corrections of the order ~ ©(m^) Eq. (H2|) coming from the right-handed coupling. 



As for MSSM— i?p[0, in addition to the charged Higgs interactions there are also 
contributions due to the sfermion, squark or slepton, mediated processes via the inter- 
actions 



(43) 
(44) 



which are obtained by expanding the additional term in the superpotential without 
-R— parity. Here the notations of f± = P±f are used for left and right-handed particles. 
It should be stressed that while both interactions in above can lead to the desired signal, 
they cannot both occur simultaneously in nature since they would lead to fast proton 
decay. 



The sfermions give a contribution like 

1 



6G, 







6G^ 



2V2G 



^ 2i3^ 3i3 p ^ i23^ »33 ^ 



(45) 



where i is the generation index with i = 1,2(1,2,3) for di{£i). For simplicity, here 
we consider only the lightest sfermions, i.e. either a down squark or slepton, in the 
calculation. With some algebraic manipulations we get the sfermion contribution as 



.MSSM-iJp 



Re 



+- 



.MSSM+i?p 



,V2GFml.F, 



1 + 




where a shorthand notation = \*\j{y2\yth) is taken. 



(46) 



4 Numerical Results 

In performing the numerical calculation, we take the central values of the known pa- 



rameters 14 



= 1.25 ± 0.15 (GeV) , = 4.1 ~ 4.4 (GeV) , mt = 173.8 ± 5.2 (GeV) , 
mw = 80.41 ± 0.10 (GeV) , Tw = 2.06 ± 0.06 (GeV) , Gp = 1-166 x 10"^ (GeV)" 

\VJ = 0.0395 ±0.0017 , sin^^iy = 0.23 . 



For X, we presume that mx sits in a region to be the order of magnitude of the top 
quark mass, and assume its decay width to be fixed to the order of Fv^sm i^i this mass 
region as well. Generally speaking the decay width would depend on the mass, but here 
for the narrow mass region in consideration, the fixed decay width assumption might 
be reasonable to certain degree. As shown in Figures |l] and ^, we see that the effect 
of varying decay width manifests significantly only in the region around threshold, i.e. 
mx ^ 110 ~ 135(GeV). 



Besides, the ratio of couplings is defined as 



= \rx\ e''^ , (47) 

since, in general, the rx could be a complex number. The angle Ox here could be an 
unremovable one relative to the phases of the 6x6 CKM matrix elements, which may 
infiuence the final result somewhat. In the case of MSSM— i?p, exploiting the bounds 



on A' and A" as given in Ref. |jT9[, we have 



< 13.61 for i = 1,2 , (48) 

r 0.06 for i = l 

< I 1.78 for i = 2 . (49) 
[ 1.32 for t = 3 

Since the lightest sfermion belongs to the third generation {i = 3), which gives the 



largest contributions in the calculation, we will take its value, = 1.32, in the whole 
sfermion mass region in consideration. To be noted that the boundarys of (^8|) and (^91) 
are obtained as taking frtj = 100 GeV and they may tend to be looser with the increase 
of mj. 

With the discussions in above on the inputs, we can now proceed the numerical eval- 
uations. Within the SM the calculation of the branching ratio is quite straightforward, 
it is 

BsMit^cT) = (y^ + 8j2ifr')'^ {^■Q±0.3±0.1)xlO-\GeVy2 

= (6.4± 1.2±0.5) X 10^^° . (50) 



Here in the evaluation, the differences between Eq. ( P^ and (|5B|), which are always 
suppressed by rhcrhr and less than 0.5%, as well among ajs, are all neglected; the i 
in the sum runs through the various color-octet states. In Eq. (^Up the errors are 
systematic uncertainties in rrit and \ Vcb\- These errors may induce about ten percent of 
uncertainties in the result as illustrated in Fig. ^. The uncertainties coming from the 
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values of mc, mr, mw and fy are apparently negligible. It is also noted that in the 
prompt T production about 30% of all may come from the feeddown of its higher excited 
states. If experiments do not separate these different sources, the present branching 
ratio should be enhanced by the same amount. 

As for the models beyond the SM, we first show the relations between the values of 
Re[Fx/ Fw^-^] and \Fx/ Fwsuf versus mx with different Tx in Fig. 0. Since the mx 
dependence exits only in these ratios, in fact Fig. ^ gives a universal feature of all X- 
mediated processes. In addition, the sizes of these two ratios are comparable and would 
decrease as mx increases, except for the region near the threshold. As a consequence of 
these facts, as long as expressions contain these ratios, the contributions are dominantly 
coming from terms including Re [Fx/FwsmI MSSM+i?p since the prefactors involving 
the ratios of couplings and /or masses are less than unity. However, for MSSM— i?p the 
term including \Fx / Fw^-^f might be dominant supposing the factors rj/\/2Gp'm^ are 
much larger than unity. 



The value of 6j, the new contribution from beyond the SM, in the MSSM+i?p is 
shown in Fig. ^ versus rrit for various mH±- For MSSM— i?p, the additional sfermion 

contributions to 5^ is displayed in Fig. ^ as a function of rrij for various 6j. Here the 

charged Higgs mass dependence of 5^ is quite weak, which is easy to figure out from 
Eq. (|Bp. 



5 Conclusions 

The top quark decay process t — > cT is investigated within and beyond the SM. In 
dealing with the heavy quarkonium production, in this paper both color-singlet and 
-octet mechanisms are employed. As a result, the decay rates are about three orders 
lower than the previous rough estimations in the SM and MSSM+_Rp. It is found 
that the decay rate even within the SM is still much larger than that of the rare top 
decays. 

As shown in Figs. ^ and ^, new interactions in models beyond the SM could enhance 
greatly the branching ratio. According to Eqs. (p7|), ( ^T]) and (^61), the branching ratio 
would be changed by a factor of (1 + i.e. B(t ^ cT) = (1 + 5j) B^uit cT), 
with 

/ 70% : 

^ [ 35000% : / ^^^^ 

in the presumed mass region given in Figs. ^ and ^, where slightly light charged Higgs 
masses are taken schematically accounting for the fact that light charged Higgs mass, 
less than 100 GeV, is still not excluded by experiments ||14[. The second enormous 
supplement of Sj stems from the enhancement of prefactor r j / ^/2G fitt-w (EH)- be- 
cause of the fact that |r^-| only runs slightly with the mj within the considered region, 
it should be pointed out that in drawing Fig. ^ we merely treat it fixed. 

Since in the next round of experiment, such as LHC, there are about 10^ ~ 10^ 
top quark events would be generated per year, it is obvious that within the SM the 
present discussed process is still at the unreachable level. Therefore, we hope this decay 
process may play a role in probing new physics in the near future. Unlike the cases 
of Quarkonium production at B-meson and Z-boson decays, here it is found that the 
color-octet contributions are negligable, although the relative coefficient enhancement 
still existing. This nature makes our calculations on a more reliable base, since it is well- 
known that large uncertainties still remain in the values of color-octet matrix elements 
among different fittings. To make a more precise estimation, the more accurate values 
on rrit and \Vcb\ are required. Last, the calculation shows that this top quark decay 
process provides a possibility of measuring the charged Higgs mass without the tan/? 
dependence, as shown in Eqs. pO|) and (p]). 
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Figure 1: \F^\/mf as a function of with = 2Twsm (solid line), = iVwsM 
(dashed line), and = Twsm (dotted line). The arrow indicates the value of m^gj^. 
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Figure 3: B{t cT) in the SM as a function of rrit with the central value \Vcb\ (solid 
hne) and the error bar coming from \ Vcb\ (short-dashed hnes). The vertical dotted lines 
indicate the experimental bounds of m^. 
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Figure 4: Sf"^ varys as a function of for mH± to be equal to 80, 110, 120, 130, 140, 
and 170 GeV, as shown in the right rectangular, respectively. The vertical dotted lines 
represent the experimental bounds of rrit. 
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Figure 5: 5j as a function of in \rj\ = 1.32 with 9j 
(dashed hne). 
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